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Abstract 



Previously published experimental work by other authors has shown 
that certain motile marine bacteria are able to track free swimming algae 
by executing a zigzag path and steering toward the algae at each turn. 
Here, we propose that the apparent steering behaviour could be a hy- 
drodynamic effect, whereby an algal cell's vorticity and strain-rate fields 
rotate a pursuing bacterial cell in the appropriate direction. Using sim- 
plified models for the bacterial and algal cells, we numerically compute 
the trajectory of a bacterial cell and demonstrate the plausibility of this 
hypothesis. 
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1 Introduction 



The experiments of Barbara and Mitchell j^] demonstrate that certain species 
of marine bacteria are able to track motile algae by a combination of reversing 
direction and steering, presumably to make use of nutrients exuded by the algae. 
Their results raise a number of interesting questions about how the bacteria are 
able to perform this feat, such as (i) how do they steer, and (ii) how do they 
decide when to reverse. In this paper we shall address the first of these questions: 
steering. We hypothesise that a bacterial cell exploits the vorticity and strain- 
rate fields generated by an algal cell in order to steer. Before presenting our 
calculations, we give a brief overview of bacterial chemotaxis. 

Bacteria play an important role in marine ecosystems. Field studies indi- 
cate that approximately half of oceanic primary production (carbon fixed by 
phytoplankton) is channelled via bacteria into the microbial loop of the pelagic 
food web (0| and references therein). This remarkable rate of bacteria-mediated 
transformation of organic matter is facilitated by bacterial motility. It is esti- 
mated that 20 to 70% of planktonic bacteria are motile [^2[ [|. Motility allows 
marine bacteria to locate sources of nutrients, for instance by colonising falling 
sinking aggregates ('marine snow') [ pTsj , or tracking motile marine algae || which 
leak organic matter p5[ . 

The most studied and well-understood motile bacterium is Escherichia coli, 
an enteric bacterium, but marine bacteria differ from E. coli in several impor- 
tant ways. First, many species of marine bacteria swim much faster than E. coli, 
with maximum speeds up to 400 /zms -1 p9| , compared to typical speeds of 
30 /jms -1 for E. coli JtJ. Second, marine bacteria are able to respond to changes 
in chemoattractant concentrations within a fraction of a second, as evidenced 
by their ability to form tight clusters f27jl , whereas E coli respond to changes 
in chemoattractant concentrations over the course of several seconds |3^] . Fi- 
nally, whereas E. coli moves in a run-and-tumblc motio n |5| , many strains of 
marine bacteria move in a run-and-reverse motion |27| . Mechanistically, 
this is explained by the fact that E. coli has several flagella which may either 
bundle together to propel the cell in a straight 'run', or fly apart to make the 
cell 'tumble' p5| , whereas many marine bacteria have only a single polar flag- 
ellum which may spin either one way or the other for a 'run' or 'reverse'. (For 
a review of other types of bacterial motions, see [^8|.) It has been suggested 
that run-and-reverse motion favoured by marine bacteria allows them to per- 
form chemotaxis (swim toward a food source) more effectively in a turbulent 
environment, compared to run-and-tumble motion p3f . Other theoretical stud- 
ies indicate that run-and-reverse may be a favourable strategy when moving up 
a uniform chemoattractant gradient in a shear flow pl[ | or when an organism is 
capable of direct gradient detection pj| . 

It seems that the run-and-reverse motion of marine bacteria is also advan- 
tageous when following a moving target such as a motile algal cell. Barbara 
and Mitchell || reported that motile marine bacteria of the species Pseudoal- 
teromonas haloplanktis and Shewanella putrefaciens were able to track motile 
marine algae of the species Pavlova lutheri by moving in zigzag paths and ap- 
parently steering, consecutively turning up to 12 times toward the algae (figure 
Q). However, P. haloplanktis and S. putrefaciens cells have a single polar flag- 
ellum [|l0| and it is not clear how they could actively steer with such an 
apparatus. Here, we suggest that the apparent steering is in fact a passive hy- 
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Figure 1: Observations of the paths of a P. lutheri algal cell being tracked by a 
P. haloplanktis bacterial cell (reproduced from figure la of 0, with permission). 
The algal cell track starts at the bottom-left of the figure. Each square represents 
a single data point; data was collected at 24 frames/sec. 



drodynamic effect, whereby an algal cell's vorticity and strain rate fields rotate 
the pursuing bacterial cell in the appropriate direction. Sections 2.1 and 2.2 



introduce simplified models for the bacterial and algal cells, based on what is 
known about them from experiments, and section 2.3 presents a scaling argu- 
ment to show that the algal cell's vorticity field has the right order of magnitude 
to account for the observed bacterial steering. Section 2.4 describes a method 
for simulating bacterial tracking, and section [3] presents simulation results and 
investigates how parameters such as swimming speeds influence the success of 
tracking. Section |] summarises the main findings and their biological relevance. 



2 Methods 

2.1 Model algal cell 

The P. lutheri algal cells used in the experiments of || had diameters of approx- 
imately 6 fim and swam at average speeds of u a i g ~ 40 /mis -1 . They changed 
direction about three times per second, and bacterial cells were able to track 
them despite these turns, but we ignore this complication and let the model 
algal cell swim in a straight line. All algae of the order Pavlovales have approx- 
imately spherical cell bodies and two flagella of unequal length, the longer of 
which is covered in fine hairs pl| . The flagella lengths of P. lutheri have not 
been measured, but in the related species P. gyrans the length of the longer 
flagellum is about 10-20 (im, and in P. mesolychnon the length of the longer 
flagellum is about 15-20 /im. The details of the flagellar beat pattern of P. 
lutheri are not well characterised, but it is known that the longer flagellum is 
directed forward during swimming. Thus, a P. lutheri cell is a 'puller', since 
it pulls itself forward with its longer flagellum, as opposed to a 'pusher' such 
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Figure 2: Simplified model of a P. lutheri algal cell, consisting of a sphere 
of radius a acted upon by a force F and a Stokeslet of strength — F located 
a distance h ahead of the sphere (the 'phantom flagellum'), such that the cell 
as a whole is force-free. In the orientation depicted here, the cell is swimming 
upward. 

as Escherichia coli whose propulsive flagella bundle lies behind the cell as it 
swims. 

In the absence of detailed knowledge of the flagellar beat pattern, we use a 
simplified model for a P. lutheri cell (figure |^). The cell body is treated as a 
sphere of radius a = 3 /im acted upon by a force F in the swimming direction 
- this force is intended to model the tension at the base of the flagellum. The 
influence of the flagellum on the fluid is modelled as a point force — F acting on 
the fluid at a distance h ahead of the cell body, with h — a equal to approximately 
half the length of the longer flagellum. (This point force 'phantom' flagellum 
is analogous to that used in [13] for modelling motile bacteria.) We choose 
h = 10 /im, but our results are qualitatively unchanged for h as low as 1 /im and 
as high as 20 /im. Note the algal cell's velocity field in this model is constant in 
time, whereas in reality it is likely to vary periodically in time with the flagellar 
beat pattern, as well as vary over a longer time-scale as the algal cell changes 
its swimming speed. 

The Reynolds number for the algal cell, based on the cell diameter, swim- 
ming speed and the viscosity and density of water, is Re rs 10~ 4 . Hence, the 
fluid flow around the cell is viscous, non-inertial, and governed by the Stokes 
equations, assuming that the fluid is Newtonian. Using the linearity of the 
Stokes equations, the velocity field around the algal cell is obtained by adding 
the solution for a translating sphere to the solution for a Stokeslet outside a 
stationary sphere p7[ ; the resulting expression is lengthy and is omitted here. 
Figure || shows the vorticity field around the model algal cell, and figure [| shows 
the direction of maximum extension rate of the strain rate field of the model 
algal cell. 
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Figure 3: Cross section of the vorticity field around the model algal cell. 
Black and white regions are where the magnitude of the vorticity field exceeds 
0.2v a i g /a, and grey regions are where the magnitude of the vorticity field is less 
than 0.2w a i g /a. Arrows show the vorticity direction. The cell body is centered 
on the origin, and the phantom flagellum is marked by a cross. Scale is in units 
of the cell radius a. 
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Figure 4: Direction of maximum extension rate of the strain rate field of the 
model algal cell. (An elongated body, such as a pursuing bacterium, will tend 
to align with the direction of maximum extension rate.) The cell body is shown 
as a black disk and the phantom flagellum is marked by a cross. Scale is in units 
of the cell radius a. In the blank regions, the direction of maximum extension 
rate is perpendicular to the page. 
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2.2 Model bacterial cell 

The two species of bacteria used in the experiments of [^| , Pseudoalteromonas 
haloplanktis and Shewanella putrefaciens, swam at average speeds of «bac ~ 
200 /ims" 1 while tracking algae. Cells from each species have a single, helical, 
polar flagellum [L4J, which may rotate either counterclockwise or clockwise 
for a 'run' or 'reverse'. Barbara and Mitchell || did not report the size of the 
bacterial cells used in their study, but cells of all Pseudoalteromonas species are 
rod shaped with width 0.2 to 1.5 /im and length 1.8 to 3 /im JnJ, while those 
of the Shewanella species are rod shaped with width 0.4 to 0.7 /zm and length 



2 to 3 /Ltm 1 36 . In each case the length of the flagellum is a few times that of 
the cell body. For the purpose of modelling we treat a swimming bacterial cell 
as a prolate spheroid of infinitesimal size which translates along its long axis 
and is rotated and advected by the velocity field of the algal cell. Ideally, one 
should incorporate the finite size of the bacterial cell and the influence of its 
velocity field on the algal cell, but to simplify the mathematics we neglect these 
complications. 

Cells of E. coli, which are a similar size to P. Haloplanktis and S. putre- 
faciens, exhibit small, random changes in direction, and these have been at- 
tributed to rotational Brownian motion (i.e. thermal collisions with molecules 
in the surrounding fluid). Berg j| estimated the the rotational diffusion co- 
efficient due to Brownian rotation to be Dr w 0.062 radians 2 s _1 for an E. 
coli cell swimming in a fluid of viscosity 2.7 cp at 32° C, and this is consis- 
tent with experimental observations ||. However, this estimate was based on 
treating the E. coli cell as a sphere of diameter 2 /im, and when one takes 
into account the stabilising effect of the flagellar bundle, the theoretically pre- 
dicted rotational diffusion coefficient is an order of magnitude smaller than the 
observed rotational diffusivity (H. Fu, personal communication), so it seems 
likely that the observed rotational diffusivity is in fact due to intrinsic 'wobbly 
swimming' rather than true (thermal) Brownian rotation (H. Berg, personal 
communication). It is therefore unclear what value of Dr is appropriate for 
P. Haloplanktis and S. putrefaciens. In section || we shall simply use a default 
value of Dr = 0.062 radians 2 s _1 , and demonstrate that the results have only a 
weak dependence on Dr anyway. 

Given our assumptions, the translational velocity of the model bacterial cell 

is 

x = v hac p + u(x,t), (1) 

where x is the position vector of the bacterial cell, Vbnc is the bacterial cell's 
swimming speed, p is a unit vector denoting the swimming direction of the 
bacterial cell, and u(x, t) is the velocity field due to the algal cell at the position 
of the algal cell. We treat Ubac as a constant, although in the experiments of 
H the bacterial cells often altered their speeds by a factor of 2 or more during 
tracking. 

With the exception of reversals and collisions (discussed later), the rate of 
change of direction of the bacterial cell is 

p = -u> x p + a p ■ E ■ (I - pp) + | x p. (2) 

The first two summands on the right hand side of (^]) represent the deterministic 
rotation of the bacterial cell by the vorticity field lj and strain-rate field E of 



7 



the algal cell p(J] , and 

a = (r, 2 l)/(r, 2 + 1), (3) 

where r\ is the ratio of the major to minor axis of the cell (slenderness ratio). 
We choose 77 = 10 as a default value (suggested as a reasonable value by PH), 
and we explore the effect of different values for r\ in section The last term 
on the right hand size of (^|) represents rotational diffusivity, with £ being a 
random, white-noise angular velocity with autocorrelation function 

&(*i)&(*a)> = AiM(*i-*2). (4) 

An important ingredient of the model is a rule to decide when the bacterial 
cell reverses direction. The details of the chemotactic response of marine bacte- 
ria are yet to be characterised p7| . For instance, it is not known whether they 
perform temporal comparisons of chemoattractant concentrations like E. coli 
|32| , whether they directly detect chemoattractant gradients, or whether they 
respond to absolute chemoattractant concentrations. Barbara and Mitchell || 
observed that during tracking a bacterial cell remained within a distance of 
about 6 /im from the algal cell. Based on this, we choose the simplest rule for 
bacterial reversals, and make the model bacterial cell reverse direction whenever 
it reaches a 'reversal distance' R « 6 /im from the centre of the model algal cell. 

One way that the bacterial cell might achieve this in reality is to reverse 
direction whenever the chemoattractant concentration drops below a certain 
threshold. The chemoattractant in the experiments is likely to have been dis- 
solved oxygen released by the algal cells (J. G. Mitchell, personal communica- 
tion), which has a diffusion coefficient of approximately 2000/im 2 s~ 1 in water 
at room temperature. Given an algal cell swimming speed of 40 /ims" 1 and an 
algal cell radius of 3/xm, the Peclet number is Pe « 0.06 <C 1. If the chemoat- 
tractant is a small organic molecule instead, then a commonly used diffusion 
coefficient is 1000/jm 2 s _1 (used for instance in |ll| |l(| |23|, ||), which gives 
Pe ss 0.12 <C 1. In either case, diffusion dominates over advection and to leading 
order the distribution of chemoattractant concentration c in the near vicinity of 
the algal cell is spherically symmetric with distribution 

c(r) = c a/r, (5) 

where cq is the concentration on the surface of the algal cell and r is distance 
from the centre of the algal cell. So, if the bacterial cell reverses direction when, 
say, r = R = 2a, the corresponding concentration threshold is cq/2. 

The final ingredient for the model bacterial cell is a rule for what happens 
when it collides with the algal cell. In reality, there will be some complicated 
interaction involving a lubrication flow between the two cells, and potentially 
also an active response from one or both of the cells. Hydrodynamic effects 
are know to cause bacteria propelled by helical flagclla to swim along curved 
paths when they are near a planar surface (|| and references therein), and 
a similar effect might occur in the interaction between the bacterial and algal 
cells. For the purpose of the model, however, we assume that when the bacterial 
cell collides with the algal cell, it simply glances off and its swimming direction 
changes so as to be parallel with the surface of the algal cell. Mathematically, 
the bacterial swimming direction pi after collision is related to the swimming 
direction p before collision by 

Pl = -7=7' ^ 
y u ■ u 



Figure 5: Idealised zigzag path of bacteria tracking algae. 



where 

u = (I - nn) -p , (7) 

I is the identity matrix, and n is the outward unit normal vector from the algal 
cell surface. 

2.3 Scaling argument 

Imagine a bacterial cell performing run-and-reverse motion in a side to side 
fashion (with respect to the algal cell's velocity) just behind a swimming algal 
cell. The algal cell's vorticity field (figure |^) rotates the bacterial cell such that 
instead of retracing the same path over and over, the bacterial cell traces out 
a curved zigzag path as depicted in figure |, with a net motion in the same 
direction as the algal cell's velocity. 

A simple scaling argument suggests that the algal cell's vorticity field alters 
the bacterial cell's path by the right order of magnitude for the bacterial cell to 
keep pace with the algal cell. If the bacterial cell is within a distance R of the 
algal cell, then the magnitude of the vorticity field it experiences is to ~ v a i g /R, 
where '~' means 'is of order'. Let W be the width of the zigzag path followed 
by the bacterial cell, and let 9 be the angle between the bacteria cell's swimming 
direction and the plane normal to the algal cell's swimming direction (see figure 
||). The angle 9 is approximately the angle through which the bacterial cell 
is rotated by the vorticity field as it swims the width of the zigzag, so 9 ~ 
wW/vh&c ~ {va\ & /vx } a. c )W/R ~ faig/wbac, since W — R. Let us assume that 
«bac S> faig, so that 9 <C 1. Then, the component of the bacterial cell's velocity 
parallel to the algal cell's velocity is UbacSinf? ~ Ubac# ~ faig- Hence the algal 
cell's vorticity field has the right order of magnitude for the bacterial cell to 
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Figure 6: Example initial condition. Circle of radius a denotes the algal cell. 
Bacterial cell initial position is random and lies a distance R from the centre 
of the algal cell. Bacterial initial swimming direction p is random and points 
inward with respect to the surface r = R (dashed circle). 



keep pace with the algal cell. Inserting numbers, w a i g rs 40/ims -1 and «bac ~ 
200 jums -1 , so 6 « 40/200 = 0.2 radians sa 11°. While this scaling argument 
focusses on the effect of vorticity, the simulation results in section || show that 
advection by the algal cell's velocity field and rotation by its strain rate field 
are also important. 



2.4 Simulation method 

It is common practice to non-dimensionalise the quantities in a problem before 
performing simulations. In the problem considered here, there are six dimen- 
sional quantities of interest (five parameters plus time): algal cell radius a, 
algal cell swimming speed v a \ s , bacterial cell swimming speed «bac, reversal dis- 
tance R, bacterial cell rotational diffusivity Z3r, and time t. In principle, by 
non-dimensionalising with respect to, say, the algal cell radius and the algal 
cell speed, one could reduce the number of quantities to just four. In practice, 
however, this is impractical, because then it is difficult to ascertain the effect 
of altering a single dimensional parameter. For this reason we shall define a 
non-dimensional time f* = iw a ig/a, but leave other quantities dimensional. One 
unit of t* then corresponds to the time it takes the algal cell to swim a distance 
equal to its own radius. 

For each simulation, the initial position of the model bacterial cell is chosen 
from a uniform, random distribution on the surface r — R, where r denotes 
distance from the centre of the algal cell. The initial swimming direction of the 
bacterial cell is chosen from a uniform, random distribution on the hemisphere 
of unit vectors that point inward with respect to the surface r — R. An example 
initial condition is shown in figure ^. The algal cell position and the bacterial 
cell position and orientation are then evolved in time, using a simple Euler time- 



step scheme, according to the rules outlined in section 2.2. The time-step At is 



chosen to be At = 0.005a/t> a i g . (Simulations were repeated using smaller time- 
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Figure 7: Example simulation trajectory of the model bacterial cell, projected 
onto a plane. The algal cell swims along the x-axis. 

steps to check convergence of results.) Brownian rotation is simulated by giving 
the bacterial cell a new swimming direction at each time-step, chosen from an 
axisymmetric distribution about the old direction, such that the new direction 
makes an angle of 2\flJ^Ki with the old direction jij. The bacterial cell is 
considered to be tracking the algal cell successfully for so long as it remains 
within the surface r — R, and is considered to have lost track of the algal 
cell if it strays outside this surface. Each simulation is terminated at either 
a pre-determined time, or when the bacterial cell loses track of the algal cell, 
whichever comes sooner. 

Unless otherwise stated, the default parameter values used in the simulations 
are: algal cell radius a = 3/xm, reversal distance R — 6/xm, algal cell swim- 
ming speed Uaig = 40/ims -1 , bacterial cell swimming speed Ubac = 200/^ms -1 , 
rotational diffusivity of bacterial cell = 0.062 /im 2 s , and bacterial cell 
slenderness ratio -q — 10. For each set of parameters, 500 simulations are run. 
We then calculate the 'empirical survival function' S(t*), defined as the frac- 
tion of simulations in which the bacterial cell is still tracking the algal cell at 
dimensionless time t*. The standard error SE(t*) of S(t*) is estimated by 



where N is the number of simulations |19[ ]. 

3 Results 

We shall first present data for the default set of parameter values listed at 
the end of the previous section, and later look at the effect of altering the 
parameter values. Figure |?] shows an example trajectory of the model bacterial 
cell, projected onto a plane, for a simulation run with the default parameter 
values. The trajectory bears similarities to the experimental results shown in 
figure 0, especially in that the primary motion of the bacterial cell is almost at 
right angles to the motion of the algal cell, so that the bacterial cell's motion 
is a zigzag. This zigzag motion is not a once-off but occurs in most of the 
simulations, as revealed by the probability density function for the angle between 
the bacterial and algal swimming directions (figure ||). For the same data set, 
figure [lO] shows the probability density of a bacterium being located at a polar 
angle 8 in a coordinate system centred on the algal cell, where 9 is measured 
from the algal cell's swimming direction. It shows that tracking most often 
occurs with the bacterial cell behind the algal cell rather than in front. 
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Figure 8: Approximate probability density function (p.d.f.) for the angle 4> 
between the bacterial cell's swimming direction p and the algal cell's swimming 
direction, for a set of 500 simulations with the default parameter values. The 
p.d.f is normalised so that in the case of isotropically distributed p it would be 
equal to unity for all (f>. Note the peak near (f> = tt/2. 
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Figure 9: Solid line shows survival function S(t*), defined as the fraction of 
simulations in which duration of bacterial tracking is at least t* , for a set of 500 
simulations with the default parameter values. Dashed lines show error bounds 
S(t*)±S E (t*). 
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Figure 10: Approximate probability density function (p.d.f.) for the polar angle 
9 describing the position of the bacterial cell relative to the algal cell, for a set 
of 500 simulations with the default parameter values. The p.d.f. is normalised 
so that in the case of an isotropic distribution it would be equal to unity for all 

e. 



Many of the tracking durations seen in the simulations are at least as long 
as those seen in experiments. Figure ^| shows the survival function and its stan- 
dard error for 500 simulations, using the default parameter values. By way of 
comparison, the P. haloplanktis and S. putrefaciens bacteria in the experiments 
of H tracked on average for 0.7 s and 1.8 s respectively, corresponding to di- 
mensionless tracking durations of t* « 9 and t* w 23. Note that the survival 
function appears to reach a steady value of S ~ 0.25, indicating that in a sig- 
nificant fraction of the simulations the bacterial cell is able to track the algal 
cell indefinitely. Extending the simulation duration up to t* = 100 reveals no 
further decay in S. Closer examination of the data reveals that in all cases 
where the tracking extends beyond t* = 30, the bacterial cell regularly collides 
with the algal cell, so re-orientations due to coll isio ns play an important role. 

The scaling argument presented in section 2JS suggests that tracking can 
occur regardless of the precise value of the algal cell's swimming speed v a i g . If 
the algal cell swims faster, for instance, then its vorticity field is proportionally 
stronger, and the consequent bacterial cell steering is more pronounced. In 
support of this, the simulation results show that the survival function is fairly 
insensitive to w a i g over a large range (20 (ims -1 < w a i g < 200 fans -1 , figure |ll| ), 
when tracking duration is measured by the dimensionless time t*. Since one 
unit of t* corresponds to the algal cell swimming a distance equal to its own 
radius, the typical distance over which tracking occurs is independent of i> a i g , 
but the absolute duration of tracking varies inversely with v a \ g . The scaling 
argument in section [2^ also suggests that tracking does not depend on the 
precise value of the bacterial swimming speed fbao an d indeed the survival 
function was found to have no statistically significant dependence on Wbac over 
the range 40/ims _1 < Vh ac < 400^ms _1 (data not shown), while holding all 
other parameters at their default values. 

Where the scaling argument fails is in its prediction that tracking is not 
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Figure 11: Survival function S(t*) for v a \ g = 20 /ims -1 (solid), w a i g = 40 /ims" 1 
(dashed), v a \ g = 80/xms _1 (dotted), and v a \ s = 200 /ims -1 (dot-dashed). Stan- 
dard errors (not shown) are of the same magnitude as those in figure ^. 
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Figure 12: Survival function S(t*) for reversal distance R = 4/im (solid), 
R = 6 /im (dashed), and R = 9 /urn (dotted). Standard errors (not shown) are 
of the same magnitude as those in figure ^. 
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Figure 13: Survival function S(t*) for default set of parameters (solid line), 
and for the same set of parameters but with rotation by vorticity and strain 
rate field turned off (dashed line). Standard errors (not shown) are of the same 
magnitude as those in figure [)| 



sensitive to the reversal distance R. As shown in figure [l2|, there is a strong 
inverse relation between 7? and tracking duration. The main reason why the 
scaling argument fails is because it does not account for advection of the bac- 
terial cell by the algal cell's velocity field, which turns out to be important. 
Taking a time-average of (|l|), the time-averaged bacterial cell velocity (x) has a 
swimming contribution i>bac(p) and an advective contribution (u), where angle 
brackets denote a time-average. For simulations with the default set of pa- 
rameters, swimming contributes ~ 30% and advection contributes ~ 70% to 
the time-averaged bacterial cell velocity. For larger R, the advection term is 
smaller, and the bacterial cell has less chance of keeping up with the algal cell. 

Given that advection contributes the dominant fraction of the time-averaged 
bacterial cell velocity, one might ask how important swimming and steering re- 
ally are. One way to test this is to artificially remove the deterministic terms 
involving u> and E on the right hand size of (||), so that there is no 'hydrody- 
namic steering'. The swimming direction p then evolves only due to collisions, 
reversals, and Brownian rotation. As shown in figure bacterial tracking 
is markedly impaired when hydrodynamic steering is turned off. Thus, even 
though advection is the dominant contribution to the net bacterial motion, on 
its own it is not sufficient to enable longer duration tracking. This concurs with 
the observations of ||, who, using latex beads as passive substitutes for bacteria, 
found that there was no tracking or entrainment of beads by any algae. 

The rate of change of bacterial swimming direction induced by the algal cell's 
strain-rate field E depends on the elongation r/ of the bacterial cell, through (||) 
and (||). A spherical cell is not rotated at all by the strain-rate field, and an 
elongated cell is rotated so as to align with the direction of maximum stretching. 
Since the direction of maximum stretching points approximately toward or away 
from the algal cell (figure ^) , one expects a more elongated bacterial cell to have 
a superior tracking performance. Figure [w] shows that this is indeed the case. 
For large rj, the precise value is not important; survival curves for different r\ 
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Figure 14: Survival function S(t*) for an elongated bacterial cell (rj = 10, solid 
line), and for an spherical bacterial cell (77 = 1, dashed line). Standard errors 
(not shown) are of the same magnitude as those in figure ^[ 



beyond rj w 5 lie almost on top of one another (data not shown) . 

As explained in section |2.2| the rotational diffusion coefficient Z?r for the 
bacterial cell is difficult to estimate, because it may be due to 'wobbly swimming' 
rather than true, thermal, Brownian rotation. Figure |l5| shows that in any 
case there is only a weak dependence on £>r. A large value of Dr does not 
significantly affect the survival function at early times t* < 20, but causes a 
slow decay in the survival function at longer times, so a bacterial cell cannot 
continue tracking indefinitely 



4 Discussion 

The simulations presented in this paper suggest that the bacterial cell steering 
seen by Barbara and Mitchell || could be a passive hydrodynamic effect. The 
model bacterial cell in the present study does not actively steer, it simply swims 
forward and back, and through its interaction with the model algal cell's velocity, 
vorticity and strain rate fields (all three of which are important) it is able to 
track the model algal cell. The simulated bacterial cell trajectories have a similar 
zigzag shape to those seen in experiments, and some of the tracking durations 
seen in the simulations are at least as long as those seen in the experiments. 
Indeed some tracking events last indefinitely (for as long as the simulations are 
run), though in reality such long tracking events would be disrupted by effects 
such as unsteadiness in the algal cell's velocity field, changes in the algal cell's 
swimming direction, or noise in the bacterial cell's chemosensory system. 

The simulations indicate that the reversal distance R (the distance from 
the algal cell at which the bacterial cell reverses direction) has a very strong 
effect on tracking duration, with smaller R enabling longer tracking duration. 
The reversal distance is an artificial construct of the model, since in reality 
the bacterial cell's chemotactic response is almost certainly more complicated 
than simply reversing direction whenever the chemoattractant concentration 
drops below a certain threshold. Even so, the inverse relation between R and 
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Figure 15: Survival function S(t*) for bacteria with rotational diffusion coeffi- 
cients Dr of radians 2 s _1 (solid), 0.062 radians 2 s _1 (dashed), 0.12 radians 2 s _1 
(dotted), and 0.2 radians 2 s^ 1 (dot-dashed). Standard errors (not shown) are of 
the same magnitude as those in figure 0. 



tracking duration shows that the ability to respond rapidly to chemoattractant 
signals and reverse direction with high frequency is crucial in order for tracking 
to occur, and this may represent an evolutionary pressure on certain marine 
bacteria toward faster response times. 

Another result of our simulations is that the dimensionless tracking duration 
(proportional to the number of body lengths the algal cell swims while being 
tracked) is largely independent of the algal and bacterial swimming speeds. This 
seems at odds with the experimental results, where || found that bacteria in- 
crease their swimming speed by roughly a factor of two when tracking (relative 
to not-tracking). If tracking occurs just as well at lower bacterial swimming 
speeds, why do the bacteria swim faster? First, the faster swimming speed in- 
creases diffusive flux, and hence nutrient acquisition, by about 30% |3|. Second, 
it could also be that faster swimming speeds result in more robust tracking when 
the algal cell changes swimming direction, a complication not included in our 
model. Third, it could be that some bacteria in the culture were committed to 
anaerobic metabolism, and only the faster-swimming aerobically metabolising 
cells performed tracking (J. G. Mitchell, personal communication). 

Perhaps the most serious limitation of the present model is that it treats 
the bacterial cell as being infinitesimally small, whereas in reality its length is 
comparable to the algal cell diameter. Thus, instead of being rotated by the 
vorticity and strain-rate fields at a point in space, in reality the bacterial cell 
'samples' these fields over its length. In defense of the present model, we point 
to the work of ]37| , who performed a combined experimental and theoretical 
investigation of the re-orientation of elongated food particles (diatoms) in cope- 
pod feeding currents. The food particles in their experiments were of similar 
length to the copepods generating the currents, but the authors still found con- 
vincing agreement with theoretical calculations that treated the food particles 
as infinitesimal. 

Given the importance of the algal cell's vorticity and strain rate fields in 
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re-orienting the bacterial cell, it is natural to ask whether ambient, turbulent 
fields might disrupt tracking. However, the strength of the strain rate tensor 
due to turbulence varies from about 1.5 s -1 in the upper mixed layer of the 
ocean under strong wind forcing to 0.005 s" 1 at the thermocline |23| ], and these 
magnitudes are small compared with the strain rate of v a \g/R w 7 s" 1 in the 
vicinity of the algal cell. Thus, turbulence is unlikely to disrupt tracking. It 
might, of course, play a role in how the bacterial cell finds the algal cell in the 
first place, which is not addressed in our model. Another consideration related 
to the strength of the algal cell's strain rate field is whether it might be strong 
enough to deform the bacterial cell's flagellum. An order of magnitude estimate 
using the bending stiffness reported by |7j] for an E. coli flagellum suggests that 
deformation could be significant. This provides another potential mechanism 
for passive hydrodynamic steering, which has not yet been explored. 

The idea that microorganisms may exploit vorticity and strain rate fields in 
their feeding behaviours is not new. In it is suggested that certain freshwa- 
ter Crustacea (cladocerans), which feed by intermittently swimming vertically 
upward and sinking down again, might enhance their prey capture rate through 
hydrodynamic effects. Their calculations show that bottom heavy motile algae 
are focussed into the wake of one of these sinking crustaceans by its vorticity 
field ('gyrotaxis'), whereupon the crustacean presumably consumes the algae 
when it swims up. In |34j] , a similar analysis is performed for the case of a bot- 
tom heavy flagellate focussed into the wake of a nonmotile, sinking algal cell. 
In the simulations of |23| , elongated model bacteria using a run-and-reverse 
('back-and-forth') motion are able to stay close to a nutrient source even at 
high ambient shear, because the strain rate field aligns the bacteria toward the 
nutrient patch. 

While the simulations presented in this paper suggest that the apparent 
steering seen by || could be a passive hydrodynamic effect, the results do not 
rule out the possibility that the bacteria are actively steering as well. The fact 
that bacteria were able to follow 'ghost tracks' some distance behind algal cells 
(mentioned at the bottom of p82 of ||) suggests that there may also be active 
steering, since vorticity and strain rate fields far behind the cell are probably 
not sufficient to rotate the cell. If P. haloplanktis and S. putrefaciens are indeed 
capable of steering, then they are not the first bacteria to be found to do so. 
Experiments indicate that another (unnamed) species of bacteria can sense an 
oxygen gradient over its body length and steer relative to the gradient in a 
continuous fashion |3^| , though it is quite different in shape from P. haloplanktis 
and S. putrefaciens. In another example, Rhodobacter sphaeroides has a similar 
shape to P. haloplanktis and S. putrefaciens, with a single flagellum, and is able 
to changes direction by altering the conformation of this helix |l[ , though there 
is no evidence of directed steering in this case. 

Finally, we should note that it has not been our intention in this report to 
make light of the work of || . Their finding of tracking bacteria is interesting in 
its own right and nutrient turnover by tracking bacteria could play a key role in 
the food web especially if, as the authors suggest, a similar mechanism applies 
to the tracking of marine snow. Rather, our intention has been to elucidate 
the intriguing possibility that not only are the marine bacteria able to perform 
tracking, but that they do so by exploiting the velocity, vorticity, and strain-rate 
fields of their quarry. 
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